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Abstract: We analyze the prospects of the detection of an 0(1) GeV neutralino dark 
matter, x?) in the Next-to-Minimal Supersymmetric Standard Model at the 14TeV LHC. 
We perform dedicated scans of the relevant parameter space of the model and find a 
large number of points where the thermal relic abundance due to such a dark matter 
is consistent with the PLANCK measurement. We note that this dark matter is highly 
singlino-dominated and is always accompanied by a pseudoscalar, Ai, with a mass around 
twice its own, which is responsible for its resonant annihilation. For two benchmark points 
from our scan, we then carry out a detector-level signal-to-background analysis of the pair 
production of a heavier higgsino neutralino and a chargino. The higgsino thus produced 
decays into the dark matter and either the Z boson or the Ai. For the Z-associated 
production of Xij we investigate the scope of the trilepton search channel. For the Ai- 
associated production mode, in order to identify the two collimated muons coming from 
the decay of the Ai, we employ an unconventional method, of clustering them together 
into one jet-like object. Using this method, for certain parameter space configurations, a 
much larger sensitivity can be obtained at the 14TeV LHC for the Aix? channel compared 
to the Zx5 channel, with an integrated luminosity of 300 fb“^. 
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1 Introduction 

The Minimal Supersymmetric Standard Model (MSSM) contains three neutral Higgs states. 
The lightest one of these, h, is one of the main candidates for the Higgs boson, Hobs; 
observed at the Large Hadron Collider (LHC) [1, 2]. The MSSM also contains a total of 
four neutralino mass eigenstates, the lightest of which, Xi^ is ™ important dark matter 
(DM) candidate, when it is also the lightest supersymmetric particle (LSP) and R-parity 
is conserved. Despite its being a theoretically sound and highly predictive model, the 
absence of any evidence of supersymmetry (SUSY) during Run-I of the LHC has brought 
the MSSM under increased scrutiny recently. This is because the only way for h to gain a 
mass close to the one measured for the Hobs is through large radiative corrections from the 
SUSY sector. The most dominant of these corrections come from the stops, which warrants 
either TeV-scale SUSY-breaking masses for them or a multi-TeV stop mixing parameter, 
At- This in turn leads to excessive fine-tuning of the model parameters such as the soft 
masses of the two Higgs doublet helds, and and the Higgs-higgsino mass paramater 
fj,, given the measured mass of the Z boson. 

The parameters of the MSSM get further constrained if one requires, in addition to h 
fulfilling the role of the observed Higgs boson, y)* to generate the observed thermal DM relic 
density of the universe, HdM) as measured by the WMAP [3] and PLANCK [4] telescopes. 
This condition is strongly dependent on the interplay between the mass and the composition 
of the XiJ as these quantities govern its annihilation rate into pairs of SM particles. The Xi 
is a pure higgsino when Mi, M 2 S> /r, with Mi and M 2 being the soft gaugino masses, and 
can only yield the correct relic abundance if /r ~ 1 TeV. If Xi is instead a bino-higgsino 
mixture, it can give sufficient Ddm for comparatively smaller masses, by annihilating via 
chargino/neutralino exchanges into pairs of vector bosons. A predominantly bino-like y)*. 
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on the other hand, can be consistent with the measured Odm for much smaller masses, 
provided one of the third-generation sfermions is not much heavier than itself. In fact, if 
such a Xi has a mass < 65 GeV, it can undergo annihilation via h in the s-channel into 
a bb pair (see, e.g., [5, 6 ]). For even smaller bino masses, the right amount of s-channel 
annihilation can occur through the Z boson. Evidently, the cross section for annihilation 
via h oi Z falls abruptly as the difference between m^o and m/i/2 or mzl‘2, respectively, 
increases. Below < 30 GeV it is thus generally impossible to obtain sufficient Odm in 

A.1 

the MSSM due to the absence of a mediator with the correct mass. 

The Next-to-Minimal Supersymmetric Standard Model (NMSSM) [7-10] contains a 
gauge singlet Higgs superfield, 5, in addition to the two doublet superhelds, Hd and Hu, 
of the MSSM. This results in two additional neutral mass eigenstates in the Higgs sector 
besides the three similar to the MSSM ones. Thus, in total there are seven Higgs bosons 
in the NMSSM; GP-even F 7 i, 2 , 3 ) GP-odd Ai ^2 and a charged pair, H^. Naturally, the two 
new Higgs bosons are dominated by the scalar component of S, which implies that their 
couplings to the fermions and gauge bosons of the Standard Model (SM) are typically 
much smaller than those of the doublet-like Higgs bosons. Their masses are thus generally 
very weakly constrained by the Higgs boson data from the Large Electron Positron (LEP) 
collider, the Tevatron or the LHG, and can be as low as a GeV or so. Similarly, the 
neutralino sector of the NMSSM also contains five mass eignestates, one more than in 
the MSSM. The fifth neutralino is dominated by the fermion component of S, commonly 
referred to as the singlino. 

In the NMSSM, the role of the Llobs can be played by either one of the Hi and 
H 2 [11-15] or even by a combination/superposition of these two when both of them lie near 
125 GeV [16-18]. Evidently, when Hohs is the Hi, only Ai can be lighter than it. Conversely, 
when ILobs is the H 2 , there exist the possibilities of either the Hi alone or both the Hi and 
Ai being lighter than it. Such light Hi or Ai can act as mediators for the annihilation of a 
Xi as light as a GeV or so in order for it to generate the correct Hdm - a scenario precluded 
in the MSSM. Note that the electroweak (EW) baryogenesis may also be accommodated 
within this framework when the Xi mass is below 10 GeV [19, 20]. Additionally, a light 
NMSSM Xi can explain the galactic centre 7 -ray excess [21, 22] observed by the Fermi 
Large Area Telesope in the presence of a light Ai [20, 23-25].^ If the Xi is even lighter, it 
has been shown [26, 27] to explain the CDMS-H event excess near 10 GeV [28]. However, 
further experimental evidence is necessary to conhrm whether any of these events are 
indeed caused by the annihilation of the DM and the jury is still out on its correct mass. 
While the direct-detection experiments such as DAMA/LIBRA [29], CoGeNT [30, 31], 
CRESST-H [32], XENONIOO [33] and LUX [34] can cover a wide range of the DM mass, it 
is extremely unlikely for them to reach below ~ 2 GeV [35] . Similarly, while the IceCube 
neutrino observatory also shows some promise for the indirect detection of the NMSSM 
DM [36], the sensitivity for a mass below 10 GeV is very poor. 

^In the case of an Hi, the s-channel annihilation is p-wave suppressed and hence generally negligible. 
Furthermore, due to the strong correlation between the masses of Hi and H 2 , as will be shown later, it is 
very difficult, although possible, to obtain an Hi as light as a GeV or so while requiring H 2 to have a mass 
near 125 GeV. 
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The Run-II of LHC can prove very crucial in this regard, owing to its potential to 
probe a very light DM, produced in association with the Z boson or the Ai in the decays 
of the heavier neutralinos. In some recent studies [37-41] it has been shown that significant 
regions of the NMSSM parameter space can be covered by the Z-associated production of 
the Xi a-t the LHC Run-H, with the Z decaying into two leptons. The X 2 3 + Xi 

channel has also shown some promise [42] even for the 8 TeV LHC, with the mass of the Ai 
lying above the production threshold. These search channels can be complemented 

by those where the light Ai is not accompanied by the DM. As long as the relic abundance 
is satisfied for a given SUSY point with a very light DM, observation of the corresponding 
Ai, which would be a likely mediator for its annihilation, could serve as a pointer. Such 
channels include, most importantly, decays of the L7obs into Ai + Ai or Ai + Z [43]. At the 
Run-I of LHC, light pseudoscalars decaying in the channel have already been probed, 

when they are produced either singly in pp collisions [44] or in pairs from the decays of a 
125 GeV Higgs boson [45]. These searches thus strongly constrain the parameter space of 
the NMSSM. Prospects for the observability of Ai at the 14 TeV LHC have been analyzed 
in detail recently in [46-48], with the measurement of the mass of Hobs providing an 
additional handle in the kinematic selection of the events. The channels investigated in all 
these analyses, however, show sufficient sensitivity only for an Ai heavier than 5 GeV. 

In this study we focus on the production of a very light, 0(1) GeV, Xi via decays of 
the heavier, higgsino-like, neutralinos of the NMSSM. We perform detailed scans of the 
NMSSM parameter space to highlight its regions where such a XiJ consistent with the 
observed relic abundance of the universe, can be obtained in the presence of a SM-like 
H 2 . We discuss some important characteristics of these parameter regions and of the Xi- 
We note, in particular, that the Xi is always accompanied by a pseudoscalar with a mass 
nearly twice its own. We then carry out comprehensive detector-level analyses of the Z- 
as well as Ai-associated productions of the Xi- For the Ai -|- Xi production channel, we 
adopt an unconventional method for reconstructing the pair of highly collinear muons that 
result from the subsequent decay of the Ai. We hnd that by using our signal reconstruction 
method the Ai-|-x5 channel can have a better signal strength at the 14 TeV LHC compared 
to the Z + Xi channel for certain specific parameter configurations. Therefore, the Ai -|- Xi 
channel can serve as a crucial probe of new physics, complementing well the Z + Xi channel 
for very a low-mass singlino-dominated DM. 

The article is organized as follows. In section 2, we will briefly discuss the NMSSM 
and its singlet sector. In section 3 we will present some details of our findings from the 
numerical scans of the NMSSM parameter space. In section 4 we will explain our signal- 
to-background analyses and in section 5 we will show their results. We will present our 
conclusions in section 6 . 

2 The singlet sector of the NMSSM 

The scale-invariant superpotential of the NMSSM (see, e.g., [49, 50] for reviews) is written 
as 

IFnmssm = MSSM Yukawa terms -|- XSHuHd + — 5*^ , (2.1) 

O 
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where A and k are dimensionless Yukawa couplings. The above superpotential observes a 
discrete Z 3 symmetry, which forbids the term present in the MSSM superpotential 

and at the same time breaks the dangerous Peccei-Quinn (PQ) symmetry [51]. Here, an 
effective /x-term is generated upon spontaneous symmetry breaking, when S develops a 
vacuum expectation value (VEV), s = (S), so that /XeS = As. The soft SUSY-breaking 
terms in the scalar Higgs sector are then given by 

b^oft = + ^^AxSHuH^ + —kAi^S^ + h.c.^ . (2.2) 

Through the minimization conditions of the complete tree-level Higgs potential, the soft 
masses mH„, ms can be traded for the respective VEVs, Vu, Vd and s, of the 

corresponding Higgs fields. 

The neutral scalar and pseudoscalar Higgs mass matrices are obtained from the Higgs 
potential evaluated at the vacuum. Diagonalization of these matrices yields the mass 
expressions for the five neutral Higgs bosons. The tree-level masses of the two lightest 
CP-even Higgs bosons, which are of our interest here, can be approximated, for moderate- 
to-large tan/3 (= Vu/vd) and EW-scale dimensionful parameters, by [52] 

~ + 4(fi;s)^ -b ksA^ 

=F \J ~ + dA^u^ [ 2 As — (H;, + ks) sin 2/3]^| , (2.3) 

where v = ~ 174 GeV. Thus there is a strong correlation between the masses of 

Hi and i/ 72 , implying that requiring one of these to lie near 125 GeV constrains the other 
also. Similarly, the tree-level mass of the singlet-like pseudoscalar can be given, assuming 
negligible singlet-doublet mixing, by the approximate expression 

9 w . . sin2/3 , , ,, 

m\ ~ A(Ha -b 4 ks)---. (2.4) 

Is 

The neutralino mass matrix in the NMSSM is written as 


/ 

All 

0 

—mw tan 9w cos 

j3 mw tan 9w sin /3 

0 \ 


0 

M 2 

mw cos /3 

—mw sin /3 

0 


-mw tan 6w cos /3 

mw cos /3 

0 

T^eff 

-XVu 


mw tan 6w sin (3 

—mw sin/3 


0 

-Xvd 

v 

0 

0 

-XVu 

-Xvd 

2ks j 


(2.5) 

with mw being the mass of the W boson and 9w being the weak mixing angle. This 
symmetric mass matrix can be diagonalized by a unitary matrix N to obtain the five 
neutralino states, Xi_ 5 , ordered by their masses, implying that the Xi is the LSP. This 
mass eigenstate is given by the linear combination 

y? = NiiB^ + Ni2W^ + /Vi3#5 + + Ni^S^. (2.6) 
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We define the gaugino fraction in x? as Zg = |A^iip + I-/V 12 P) the higgsino fraction as 
Zh = I-/V 13 P + I-/V 14 P and the singlino fraction as Zg = iWsp. The composition of x? 
thus depends on the values of the various model parameters appearing in the above mass 
matrix. 

In the MSSM, the fifth row and column, corresponding to the singlino and hence 
containing the A-dependent terms, do not exist. There, the smaller the \i compared to the 
min[Mi,M 2 ], the larger the higgsino fraction in Xi- Since the mass of the higgsino-like 
chargino, which is the lighter Xi^ in this case, is also proportional to fi, the value of this 
parameter is bounded from below to around 100 GeV by the non-observation of a chargino 
at the large electron positron (LEP) collider. On the other hand, in order to avoid excessive 
fine-tuning of the model parameters for obtaining the correct Higgs boson mass, /i ought 
not to be too large. However, a purely higgsino LSP does not give the correct thermal 
relic density for masses below ^ 1 TeV, as the annihilation cross section becomes too large. 
Therefore, the only way for the x? in the MSSM to give good relic abundance with a mass 
0(10) GeV is to have a sizable bino component [ 6 , 53]. 

In the NMSSM, the presence of the singlino leads to some additional possibilities in 
the context of DM phenomenology. In the limit ^eS ^ min[Mi,M 2 ], the term [Ad^olss = 
2ks = in eq. (2.5) implies that the LSP is singlino-dominated for 2k/X < 1. The 

singlino fraction in ^ can be increased by reducing k. One could alternatively increase A, 
while keeping //efr fixed, to reduce the size of the [Adpjss term, but this also enlarges the 
sizes of the off-diagonal mixing terms. Eq. (2.4) shows that the mass of Ai also scales with 
Ks. Thus a light singlino-like Xi can be naturally accompanied by a light Ai. Note that the 
strong correlation, eq. (2.3), between the H 2 mass, which is required to be around 125 GeV, 
and the Hi mass generally prevents the latter from taking very low values, although it also 
scales mainly with ks. Importantly though, even if Hi has a mass close to 2x5, A? 
annihilation via s-channel Hi is p-wave suppressed, which would make it extremely difficult 
to generate the correct thermal relic abundance. We therefore focus only on a light Ai 
here, which can result in the singlino-like x 5 yielding the correct relic density even with its 
mass around a GeV. 

As noted in the Introduction, the light x? can be produced in the decays of the heavier 
neutralinos. In particular, the X 2 ^cid X 3 as well as the xt in this scenario with a light 
singlino-like DM are predominantly higgsinos [42] . Their main production channel is pp —>• 
X 23 + xfj which is followed by the decays X 23 ZjAi -|- x5 -|- tiT nnd xf 

+ Xi where implies missing transverse energy. The complete processes 

are shown in Fig. 1, (a) and (b). The .^-|-Xi decay channel of the X 2 3 is by far the dominant 
one, while X 2 3 + Xi is generally suppressed. The Z + W + $t production, which 

results in trileptonic final states, is therefore the preferred search channel for the higgsinos, 
also since it has minimal dependence on m^o. Dedicated searches carried out by both 
the CMS and ATLAS collaborations [54-56] have either already put strong constraints on 
significant regions of the NMSSM parameter space or are likely to cover large portions of 
it, where the x? can be as light as a GeV, at the LHC Run-H. 

It is, however, important to note that, for sufficiently large values of the parameter A, 
the branching ratio (BR) for the X 2 3 ^i+X? decay channel, while still subdominant, can 
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Figure 1: Production process for the X 2 , 3 xf pair, with the X 2,3 decaying via (a) Z + Xi 
and (b) + x?- 


become sizable. For Ai masses of interest here, i.e., below the threshold, the decay 

mode with the largest BR is Ai —)• cc. The main leptonic decay channel, Ai —)• is 

highly subdominant, with its BR never exceeding 9%, as we shall see below. This value 
is still larger than twice that of the BR(Z —)• (~ 0.034). But one challenge that 

arises in this context is the reconstruction of the very light Ai from two highly collinear 
muons it decays into and the isolation of this signal from the background. As long as this 
complication can be overcome, the two observations above imply that the Ai + y® search 
channel can possibly complement the Zxi channel for certain specific NMSSM parameter 
space points. In the following sections we will analyze such parameter combinations and 
also introduce our method for reconstructing the very light Ai from two collinear muons. 
We should point out here that while in principle the same method can alternatively be 
used to probe these Ai’s in the decays of the heavy CP-even Higgses, our requirement of 
the presence of in the final state, as an indication of light DM, renders these channels 
irrelevant here. 

3 Parameter space scan and constraints 

Due to the presence of the additional singlet superfield, the NMSSM contains several new 
parameters besides the 115 or so of the MSSM at the EW scale. However, assuming the 
sfermion mass matrices and the scalar trilinear coupling matrices to be diagonal reduces 
the parameter space of the model considerably. In this study, since our main focus is the 
higgsino-singlino sector, we further assume the following universality conditions. 



(3.1) 
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Parameter 

Scanned range 

Mo (GeV) 

500 - 2000 

Mi /2 (GeV) 

300 - 1000 

Ao (GeV) 

0 - 4000 

9es (GeV) 

100 - 300 

tan /3 

6-25 

A 

0.01 - 0.4 

K 

10 "® - 10 -^ 

AIa (GeV) 

0 - 5000 

A« (GeV) 

-100 - 0 


Table 1: Ranges of the NMSSM input parameters scanned for obtaining Xi with a mass 
below 2 GeV. 


where Mq^ 2 , 3 ^ -^[/i, 2.3 5 -^£> 1 , 2 , 3 ) -^^ 1 , 2.3 -^^ 1 , 2,3 the soft masses of the sfermions, 

-^^ 1 , 2,3 those of the gauginos and At^^r the trilinear Higgs-sfermion couplings. Along with 
Mq, Mi /2 Aq, the model then contains A, k, /TeS, tan/3, Ax and as the only free 
parameters, which are input at the SUSY-breaking scale, Msusy = > with ^ 

being the physical masses of the two stops. 

We performed a scan of the above mentioned nine parameters using the nested sampling 
package MultiNest-v2.18 [57], which was interfaced with NMSSMTools-v4.5.0 [58-61] 
for calculation of the SUSY mass spectrum and BRs for each model point sampled. The 
scanned ranges of the parameters are given in table 1 and were chosen, based loosely on 
the Endings of [42] , to maximally yield points with m^o < 2 GeV as well as the H 2 having 
a mass near the one measured for the iifobs at the LHG. The H 2 was additionally required 
to have SM-like signal rates. The signal rate, obtained from NMSSMTools as an output, 
is defined, for a given decay channel X, as 


— ^{99 H 2 ) X BR(ff2 —^ X) 

^{99 ^sm) X BR(/ism — X) 

where hsu is the SM Higgs boson with a mass equal to • 

A very light Ai is subject to constraints from direct collider searches as well as from 
flavor physics [62]. The program NMSSMTools intrinsically takes into account the ex¬ 
clusion limits on pseudoscalars from LEP and BaBar as well as from the 4/r searches at the 
CMS [45]. Additionally, It ensures that the LEP limit on the xt mass and the perturba- 
tivity constraints on the various Higgs boson couplings are satisfied and that the x? is the 
LSP. We further required each scanned point to satisfy the following constraints: 


• 2.63 X lO-"^ < BR (B ^ Xsy) < 4.23 x lO”'^, 

• 0.71 X 10"^ < BR(Bu ^ Ti/) < 2.57 x 10"^, 

• 1.3 X 10-9 < BR(Bs ^ < 4.5 x lO'^, 


- 7 - 







• 0.107 < < 0.131, 

• 122 GeV < mH 2 < 128 GeV. 

The 6 -physics observables above were calculated using the package SuperIso-v3.4 [63] 
and the relic density using the MicrOMEGAs-v4.1.5 [64] package. The allowed range of 
assumes a ±10% theoretical error around the central value of 0.119 measured by the 
PLANCK telescope [4]. Similarly an error of ±3 GeV is allowed in the theoretical estimation 
of the H 2 mass, given the Lfobs mass measurement of 125 GeV at the LHC. Consistency 
with the LEP and LHC exclusion limits on the non-SM-like Higgs bosons, including the 
Ai, was further ensured by testing each point that passed the above constraints with the 
program HiGGSBouNDS-v4.2.0 [65-68]. 

In Fig. 2(a) we show the ranges of the Ai and Xi masses obtained for the good points 
from our scan. The heat map corresponds to the Xi relic density. One can see a strong 
correlation between m^o and ruAi, with the former almost always being half of the latter 
in order for the resonant annihilation of the Xi via the Ai. The tiny Xi masses are a result 
of vanishing k as seen in Fig. 2(b), implying an almost PQ-symmetric model. The heat 
map in the figure shows the distribution of the singlino fraction in the Xi which increases 
as K decreases and is always larger than 0.75. A large singlino component implies a small 
higgsino fraction, which is necessary to prevent an under abundance of the from too 
much annihilation. Fig. 3(a) illnstrates that the Ai is restricted to a lower mass, needed for 
satisfying the relic density constraint, by adjusting to a narrow range of low negative 
values. This is in agreement with eq. (2.4), along with the fact that A, illustrated by the 
heat map in the fignre, also generally tends to be small. In Fig. 3(b) the distributions of the 
parameters Ax and //eff are shown against the Hi mass range obtained in the scan. We see 
that when both Ax and tan/3 are small rriHi is low, while its maximum value, ~ 45GeV, 
is obtained for the largest allowed values of Ax, with tan/3 > 10. 

Fig. 4(a) shows that the mass of the X 2 almost equal to the parameter /Teff as long 
the parameter M 1 / 2 , given by the color map, approaches its maximum allowed value. As 
the splitting between M 1/2 and /Teff decreases, the gaugino-higgsino mixing increases, which 
results in somewhat lowering m^o relative to /Teff- Fig. 4(b) similarly shows that, since the 
X 3 and the xf S'l's higgsino-like as well, their masses are also proportional to //eff- Moreover, 
the xf is always almost mass-degenerate with x®, while X 3 is generally heavier than both 
of them, but only slightly so.^ Fig. 5(a) illustrates that the BR(x 2 ^ 1 X 1 ) increases with 
A, shown by the heat map. This is because, for vanishing singlet-doublet mixing in the 
pseudoscalar mass matrix as well as singlino-higgsino mixing in the neutralino mass matrix, 
the coupling between the Aj, Xi and X 2 ) given by eq. (A.14) of [49], can be approximated 
by the simple relation 

^^ 1 X 1 X 2 ~ ~ 2-\/2k) . (3.3) 

^Due to the nature of the neutralino mass matrix, given in eq. (2.5), both positive and negative values 
of the X 3 mass are possible, although we show here only negative valued solutions. Positivity of m~o can 
be ensured a priori by a phase transformation. 
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(a) (b) 

Figure 2: (a) The distribution of the Xi mass vs. that of the Ai mass, with the heat 
map corresponding to the relic abundance, (b) The dependence of m~o on the parameter 
K. The heat map shows the singlino fraction in the Xi- 
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Figure 3: (a) The mass of Ai as a function of the parameters 4^ and A. (b) The mass of 
Hi vs. the distributions of the parameters Ax and tan/3. 


Still, the BR(x 2 ^iX?) exceeds 10% only for a few points, reaching as high as about 30% 
for a couple of them. On the other hand, the BR(x 2 given on the x-axis, drops 

quite sharply with increasing A and can reach as low as about 30%. This is due mainly 
to the fact that the BR(x 2 ^ -f^2Xi) (riot shown here) also rises, much more abruptly 
than the BR(x 2 ^iXi); as A increases. Fig. 5(b) shows that the BR(x 3 —>■ Ai^) never 
exceeds 9% while the BR(a 3 —> .^Xi) never falls below 50%, implying a relatively weaker 
dependence on A. The heat map in the figure shows that the maximum BR(4i —> 
achievable is 9%. Note that the BR(xf —^ H^^Xi) is unity for all the points shown in these 
figures. 

We point out here that we also carried out a test scan of the NMSSM with partial 
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Figure 4: (a) The X 2 mass as a function of the parameter figg. The color map shows the 
dependence on Mi/ 2 - (b) The distribution of the X 3 and xf masses vs. ^gg, shown by the 
heat map. 



(a) (b) 

Figure 5: (a) The BR(x 2 ^xT) ^be BR(x 2 ^ 1 X 1 ) as functions of the parameter 
A. (b) The BR(x 3 —^ ^Xi) vs. the BR(x 3 —>• Aix^), with the heat map corresponding to 
the BR(Ai —)• 


universality at the grand unihcation (GUT) scale. In this model the unified scalar mass 
parameter, mo, the universal gaugino mass, mii 2 , the universal Higgs-sfermion trilinear 
conplings, oq, as well as the paramaters a\ and are input at the GUT scale. The EW 
scale values of the individual soft scalar and gaugino masses and of all the Higgs trilinear 
couplings are obtained by the running of these parameters using the renormalization group 
equations. The other Higgs sector parameters. A, k, figg and tan/3 are input at Msusy- 
But, owing particularly to the GUT scale definition of the parameters ax and a^, severe 
fine-tuning is necessary in this model for the scan to find the right combinations of these 
parameters at MgusY that yield both H 2 and Ai with the desired masses. 
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Selection 

SRnoZa 

SRnoZb 

SRnoZc 

SRZa 

SRZb 

SRZc 

tosfos 

< 60 

60-81.2 

< 81.2 or > 101.2 

81.2-101.2 

81.2-101.2 

81.2- 101.2 

tr 

> 50 

> 75 

> 75 

75 - 120 

75 - 120 

> 120 

]\4j' 

- 

- 

> 110 

< no 

> no 

> no 

PTi^s) 

> 10 

> 10 

> 30 

> 10 

> 10 

> 10 

SR veto 

SRnoZc 

SRnoZc 

- 

- 

- 

- 


Table 2: Selection requirements for the six signal regions defined for the trilepton searches 
by the ATLAS collaboration. All the dimensionful parameters in rows 2 — 5 are in units 
of GeV. 


As noted above in Fig. 3(a), at Msusy is restricted to a narrow range of values, 
which would imply an even smaller set of its possible values at the GUT scale. Furthermore, 
the SUSY-preserving parameter k is also very small due to our requirement of ~ 
1.5 GeV, though the approximate PQ symmetry this results in still avoids the cosmological 
constraints on the PQ axion [69]. Keeping /Xeff around the EW scale, k is typically smaller 
than 0.01 A. Despite all these limiting conditions, some points with the right masses of the 
Xi and Ai were indeed found. However, they did not cover, for example, the wide range 
of the BR(x 2 —)■ AiXi) ™ 5(a), which stretches beyond 0.1. We therefore opted 

for the 9-parameter EW-scale NMSSM for this study. 

4 Very light DM via higgsino decays at the LHC 

Eor our signal-to-background analyses, we selected two benchmark points, BPl and BP2, 
out of the good points from the NMSSM parameter space scan. BPl is chosen such that 
the BR(y 2 3 -^Xi) is sufficiently large, while BP2 has a relatively enhanced BR(y 2 3 —>■ 

AiXi)- The parton-level signal and background events for these points were generated with 
MADGRAPH_aMG@NLO [70] and passed to Pythia 6.4 [71] for hadronization. 

4.1 The trilepton channel 

As noted earlier, the ATLAS and GMS experiments have separately performed searches 
for trileptons (Si) [54-56] resulting from the X 2 sxf pah production. In the ATLAS 3i 
search [54], which is the one we will consider in the following, six signal regions (SRs) are 
defined in terms of the invariant mass of two same flavor leptons with opposite sign (SEOS), 
"iSFOS- These regions also depend on the momentum, ^^(£ 3 )) of third lepton, ^ 3 , that is 
left after requiring two SFOS leptons to reconstruct msFOS) $t th® transverse mass, 
Mt = ^ 2 .^ 2 ’Pt( 4 ) (1 ~ cos where is the azimuthal angle between 

and the ^ 3 . The signal regions are divided into three ‘Z-enriched’ ones, SRZ{a,b,c}, where 
msFOS lies within 10 GeV of mz-, and three ‘Z-depleted’ ones, SRnoZ{a,b,c}, where mspos 
lies outside this mass window. Table 2 further shows the selection requirements for each 
of these six regions. 

The irreducible backgrounds include di-boson, tri-boson and tiWjZ productions, all 
of which can have three or more leptons and in the final states. The ZZ and ZW^ 
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backgrounds are by far dominant over the tiW/Z one [54], Among the reducible back¬ 
grounds are included top quarks produced singly or in pairs, WW and W ov Z bosons 
produced in association with jets or photons. Among these the ti background is highly 
dominant. For each BP, the cross section for the signal process was calculated at next- 
to-leading order (NLO) using Prospino-v2.1 [72]. We then first generated the event files 
corresponding to y/s = 8 TeV for the signal process and passed these to the public package 
CheckMATE-v1.2.0 [73] for testing against the current LHC limits from the trilepton 
searches. In CheckMATE the signal regions given in table 2 have been defined and the 
corresponding backgrounds from the ATLAS experiment implemented. For testing a model 
point it therefore calculates the signal efficiency for each region, after ATLAS detector sim¬ 
ulation with DELPHES 3 [74]. 

After confirming that the given BP is not excluded by the available data, we proceeded 
to the future 3£ search at the 14 TeV LHC. We generated the signal event files for y/s = 
14 TeV and passed these to CheckMATE again. In this way we obtained the signal 
efficiencies, after multiplying the NLO cross section with an assumed integrated luminosity, 
£, of 300fb“^ (be., design luminosity of the LHC) to get the number of signal events in 
each signal region. 

As for the backgrounds, we only simulated the three dominant ones, ZZ, ZW^ and tt, 
for the 14 TeV LHC. We used MadGraph to generate the background events and passed 
these to CheckMATE to get the cut efficiencies for all the backgrounds in each signal 
region. After multiplying the NLO cross section [75, 76] and the luminosity, we got the 
number of background events in each signal region. 

4.2 Collimated muons from an Ai 

Due to the smallness of the Ai mass of our interest here, the muons it decays into are 
highly collinear. In order to isolate such muons, we employ the technique of clustering 
them together into one object, ^co\- This method, similar in concept to the construction 
of a ‘lepton-jet’ [77], has recently been shown in [78] to be very effective for probing highly 
mass-degenerate higgsinos. 

For using this method, the signal events generated for BPl and BP2 were passed to 
Pythia 6.4 for hadronization and subsequently to DELPHES 3 for jet-clustering via the 
anti-/cT [79] algorithm using FastJet-v3.0.6 [80]. The object ^coi is constructed as follows. 

1. Require the transverse momentum, pT, larger than lOGeV for each muon in the 
signal. In addition, impose the cut < 5 GeV on the invariant mass of the muon 
pair. 

2. Define /gum as the scalar sum of the transverse momenta of all additional charged 
tracks, each with px > 0.5 GeV, within a cone centered along the momentum vector 
of pcoi and satisfying Ai? = 0.4. Impose igum < 3 GeV. 

The main backgrounds, containing two collinear muons along with a third lepton and 
$Ti include 1V(—)• and Z(—wherein the ^col comes from the photon, 

and Whh, when one of the 6 -jets produces the pcoi- The Z 7 * background fakes the signal 
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process when one of the two leptons escape undetected. Note that while the tt background 
mentioned above for the 3i search is also relevant for this signal process, it becomes negli¬ 
gible here owing to the requirement of the two hnal state muons being highly collinear. In 
order to isolate these backgrounds, we implement the following cuts. 

1. Since the Ai resulting from the higgsino decay is highly boosted, we expect the ficoi 
from its subsequent decay to have a large pT- We therefore require PT(/ieoi) ^ GeV. 
We also require Pt((.^) > 20 GeV. 

2. In order to reduce the background a large is required, so we add the cut > 
50 GeV. 

3. In the IV 7 * background the Mt distribution has an end point around the W boson 
mass, which leads us to impose Mt > 80 GeV. 

4. Our signal would appear as a narrow peak in the distribution. Hence, we 

impose a narrow cut width, 5 x ta^-, where cjai = 0.26 + O.OOlSruAi, around niAx- 
This parametrization of aA^ follows the mass resolution study of the J/\k resonance 
in [81]. We also remove the J/'k resonance region (3.0GeV < < 3.2GeV). 

To get a sufficient number of Monte Carlo events in the kinematic regime of our interest, 
we require > 1.5 GeV and PT(ii) > 10 GeV at the parton level for the IT 7 * and Z 7 * 
backgrounds. For the Whh background, we additionally require the pT of the 6 -jet to be 
larger than 30 GeV. 

5 Results and discussion 

In table 3 are recorded some specifics of the two BPs used for the signal-to-background 
analyses. The consistency of each of the H 2 signal rates given in the last three rows of 
the table is to be checked against the experimental quantity ^^^h 

corresponding channel X. Note that this comparison assumes that the inclusive pp cross 
section for H 2 production can be approximated by the dominant gluon-fusion production 
cross section. Note also that, since the WW and ZZ decays of H 2 are proportional to the 
same coupling, NMSSMTools provides a unique value of the signal rates for these two 
channels, which we denote by in the table. The most recent measurements of by 
the GMS (ATLAS) collaboration(s) read [82-85] 

p^^ = 1.13 ±0.24 (1.17 ±0.27), 
pZZiww) ^ 1,0 ±0.29 (1.09l°;23), ( 5 .I) 

= 0.91 ±0.28 (1.4l[]j). 

Using each of the two analysis methods described in the previous section we calculated 
the number of signal events, S, and that of background events, B, for each BP at the LHC 
with y/s = 14TeV and L = 300 fb“^. In table 4 we provide the S and B corresponding 
to each of the six signal regions in the ATLAS 3i search. The total signal cross sections 
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BPl 

BP2 

Model parameters 

Mo (GeV) 

1951.1 

1826.0 

Mi /2 (GeV) 

892.24 

929.2 

Zlo (GeV) 

2462.2 

2626.2 

ties (GeV) 

191.34 

164.52 

tan/3 

14.056 

19.785 

A 

0.0814 

0.3102 

K, 

0.0002 

0.0008 

Aa (GeV) 

4080.2 

3596.7 

41. (GeV) 

-3.6681 

-6.8466 

Masses 

m^o (GeV) 

1.0025 

1.4081 

(GeV) 

189.09 

170.13 

m^o (GeV) 

-201.67 

-182.27 

m~± (GeV) 

194.97 

167.72 

(GeV) 

2.1776 

2.9856 

mH 2 (GeV) 

124.12 

125.79 

Branching Ratios 

BR{xl ^ Zx\) 

0.634 

0.603 

BR{xl ^ Aix?) 

0.004 

0.089 

BR{xl ^ Zx?) 

0.736 

0.704 

BR{xl ^ Aix?) 

0.004 

0.081 

BR{Ai —>• 

0.039 

0.087 

H 2 signal rates 


0.998 

0.901 

R^^ 

0.996 

0.885 


1.003 

0.847 


Table 3: Properties of the two benchmark points used for the signal-to-background anal¬ 
yses. 

obtained are 24.3 fb and 3.93 fb for BPl and BP2, respectively. In the Z-enriched region, 
ZW^ production dominates the total background. In the Z-depleted region, a comparable 
contribution is obtained from the tt background. One can notice in the table that for both 
the BPs, the highest sensitivity is obtained in the signal region SRZc. 

In the pcoi channel, the cross sections for the signal and background processes are given 
in table 5 for the BPl, along with the efficiency of the cuts and the effective cross sections 
and numbers of events after implementing these cuts. One notices that the Wbb is by far 
the largest background. It, however, gets greatly reduced by the cuts, after which IP 7 * 
takes over as the most dominant background. The corresponding quantities for the BP2 
are listed in table 6 . Due to the different masses of the Ai and Xi obtained for the two 
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Background or signal 

SRnoZa 

SRnoZb 

SRnoZc 

SRZa 

SRZb 

SRZc 

ZZ events 

410 

59 

10 

280 

39 

12 

ZW^ events 

1391 

595 

71 

6850 

661 

189 

tt events 

1715 

401 

62 

272 

178 

19 

All background events 

3516 

1055 

143 

7402 

878 

220 

BPl signal events 

12 

37 

19 

191 

134 

130 

BP2 signal events 

20 

46 

18 

270 

144 

96 


Table 4: The number of background and signal events at the 14 TeV LHC for C = 300 fb ^ 
in each of the signal regions of the ATLAS 3£ search. 



BPl 

Wy* 

Zy* 

Wbb 

Cross section (fb) 

0.178 

246.9 

10.0 

3770.0 

Cut efficiency 

0.123 

2.15 X 10"^ 

6 X lO"*^ 

1 X 10"® 

Effective cross section (fb) 

0.022 

0.053 

0.0006 

0.003 

No. of events 

6.6 

15.9 

0.18 

0.9 


Table 5: The backgrounds and the signal for the BPl in the //col search channel at the 
14 TeV LHC for C = 300 fb'^ 


BPs, different sets of cuts need to be implemented for each of them. 

We quantify the strengths of the two analyses in terms ol S/B for comparing them 
against each other. This quantity is listed in table 7 for the two BPs in each of the search 
channel. The 3^ analysis gives a slightly higher S/B for the BPl compared to that for the 
BP2. On the other hand, S/B for the BP2 in the ^coi analysis is considerably larger than 
the S/B obtained for the BPl in each of the two analyses. This is due, evidently, to the 
sizable BR(x 2 —^ ^iXi) BR(Ai —)■ as noted in table 3. In addition, the cut 

efficiency for the signal is much higher while that for the Wj* background is much lower 
in the case of the BPl. 

For a more realistic analysis of the prospects of a signal process though, the statistical 
and systematic uncertainties in it also need to be taken into account. Hence, for each 
BP we have also provided in table 7 the statistical significance, given by the approximate 
formula. 




S 

VB + (eB)2 ’ 


(5.2) 


where the systematic uncertainty is given by the fraction £ of the background. From the 
ATLAS 3i search [54], we note that the systematic uncertainty is 21% for the SRZc signal 
region, where the highest sensitivity is achieved, as seen above. We expect this number 
not to vary considerably at the 14 TeV LHC and hence use e = 0.21 in our estimation 
of Z for the 3i channel. For the W + ficoi + channel, since there is no experimental 
analysis available, the systematic uncertainty has to be estimated roughly. There are two 
major sources of this uncertainty: the reconstruction of the ficoi, in which case it is around 
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BP2 

IU 7 * 

Z 7 * 

Wbb 

Cross section (fb) 

3.93 

246.9 

10.0 

3770.0 

Cut efficiency 

0.050 

5.3 X 10"^ 

3 X 10-5 

1 X 10-6 

Effective cross section (fb) 

0.197 

0.013 

0.0003 

0.003 

No. of events 

59.1 

3.9 

0.09 

0.9 


Table 6: The backgrounds and the signal for the BP2 in the /Xcoi search channel at the 
14TeV LHC for C = 300fb"!. 


Point 

S/B in analysis 

Z (a) in analysis 


3£ (SRZc region) 

/^col 

3i (SRZc region) 

l^col 

BPl 

0.591 

0.42 

2.7 

1.2 

BP2 

0.436 

15 

2.0 

27 


Table 7 : Results in the two analyses methods for the benchmark points. 


5% [45], and that of the £ 3 , where it is less than 2 % [ 86 ]. As a conservative estimate, which 
also allows a direct comparison between the 3£ channel and this channel, we set £ = 0.21 
here also. This results in = 27cr in this channel for the BP2, as seen in table 7, which 
is much higher than the estimated Z for the same point in the 3^ channel. 

There is, however, a caveat here. As noted from table 6 , B for the BP2 is much 
smaller than S, resulting in a large S/B. In such a case, the approximate expression for Z, 
which assumes S' ^ B, is in principle not valid [87]. While, for a consistent treatment of 
the systematic uncertainties between the two search channels, we retain this approximate 
expression for the /icoi channel also, we emphasize that the given S/B values be considered 
a much more accurate estimate of the strength of each channel at the 14TeV LHC. 

6 Conclusions 

In this article we have discussed an 0(1) GeV neutralino DM in the NMSSM and its 
detectability at the LHC. Despite being very light, such a singlino-like DM can generate 
thermal relic abundance of the universe consistent with the PLANCK measurement, owing 
to the existence also of a singlet-like pseudoscalar, Ai, with a mass around twice the DM 
mass. A very light DM giving the correct relic abundance is impossible to obtain in the 
MSSM, and thus its detection will provide a clear indication of physics beyond minimal 
SUSY. We have noted that the current direct and indirect detection facilities have very 
poor detection prospects for such a light DM with and hence its dedicated searches at 
the LHC can prove very crucial. We have therefore studied in depth the prospects for its 
observation at the 14TeV LHC. 

By performing a through scan of the NMSSM parameter space, chosen taking into 
account the analytical structure of the neutralino mass matrix, we found a variety of its 
configurations where an 0(1) GeV DM can be obtained. We then carried out detector-level 
analyses of two of the main production modes of such a DM. In both these modes, the DM 
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is produced in the decays of a higgsino-like heavier neutralino, X2/3) which itself is produced 
in a pair with the lightest chargino. The X 2/3 then decays into either Z + DM or Ai +DM. 
The former channel, combined with the chargino decay, results in a trilepton + final 
state, for which dedicated searches are already being performed by the CMS and ATLAS 
collaborations. In the latter channel the final state comprises of a pair of muons and a 
third lepton along with the DM. 

In the X2 3 h)M + Ai channel, the two muons that the very light Ai decays into 
are highly collinear. Therefore, this channel can not be probed using the usual muon 
identihcation criteria. For this reason, we have adopted the technique of grouping the two 
muons together into a single jet-like object by applying certain unconventional rigorous 
cuts. By implementing this method on two benchmark points from our scan, we have 
found that this channel can have a signal strength comparable to that of the trilepton 
channel at the 14TeV LHC with L = 300 fb~^. In fact, for one of the two points, wherein 
the BR(x 2 3 DM-I-Ai) is significant and the BR(Ai —>• is maximal, the obtained 

SjB for this cannel is much larger than that for the trilepton channel. We thus emphasize 
that dedicated searches in this channel may prove very crucial for the discovery of a very 
light SUSY DM within a few years of the current LHC run. 
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